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Motivated by the string corrections on the gravity and electrodynamics sides, we consider a 
quadratic Maxwell invariant term as a correction of the Maxwell Lagrangian to obtain exact solutions 
of higher dimensional topological black holes in Gauss-Bonnet gravity. We first investigate the 
asymptotically flat solutions and obtain conserved and thermodynamic quantities which satisfy 
the first law of thermodynamics. We also analyze thermodynamic stability of the solutions by 
calculating the heat capacity and the Hessian matrix. Then, we focus on horizon-fiat solutions with 
adS asymptote and produce a rotating spacetime with a suitable transformation. In addition, we 
calculate the conserved and thermodynamic quantities for asymptotically adS black branes which 
satisfy the first law of thermodynamics. Finally, we perform thermodynamic instability criterion to 
investigate the effects of nonlinear electrodynamics in canonical and grand canonical ensembles. 


I. INTRODUCTION 


Nonlinear theories have been studied extensively in the context of various physical phenomena. In other words, 
because of the nonlinear nature of most physical systems, linear theories could not precisely describe the experimental 
consequences and we must inevitably consider the nonlinear models. 

Classical electrodynamics, whose field equations originate from the Maxwell Lagrangian, contains various problems 
which motivate one to consider nonlinear electrodynamics (NLED). One of the main problematic items of Maxwell 
theory is infinite self-energy of a point-like charge. The first successful Lorentz invariance theory for solving this 
problem was introduced by Born and Infeld [l[. After Born-Infeld (BI) theory, although some different NLED have 
been introduced with various motivations, their weak field limits lead to Maxwell theory Amongst the NLED 

theories, the so-called Bl-types, whose first nonlinear correction is a quadratic function of Maxwell invariant, are 
completely special In addition to the interesting properties of BI-tw)e NLED theories , it may be shown that 

these theories may arise as a low energy limit of heterotic string theory Q , which leads to increased interest. Taking 
into account the first leading correction term of Bl-type theories in which originated from the string theory, one can 
investigate the effect of NLED versus Maxwell theory [^, [l^ . 

In this paper, we consider Gauss-Bonnet (GB) gravity as a natural generalization of Einstein gravity with a quadratic 
power of Maxwell invariant in addition to the Maxwell Lagrangian as a source. On the gravity point of view, string 
theories in their low energy limit give rise to effective models of gravity in higher dimensions which involve higher 
curvature terms, while from the electrodynamics viewpoint, a quadratic power of Maxwell invariant supplements to 
the Maxwell Lagrangian with the development of superstring theory. In the weak field limit, GB gravity and the 
Lagrangian of the mentioned NLED reduce to Einstein gravity and Maxwell Lagrangian, respectively. 

On the other hand, it has been shown that there is a close relationship between black hole thermodynamics and 
the nature of quantum gravity. For black holes to be physical objects, it is necessary for them to be stable under 
external perturbation. There are two kinds of stability: physical (dynamical) stability and thermodynamical stability, 
which is of interest in this paper. The most interesting parts of black hole thermodynamics are thermal instability 
and Hawking-Page phase transition ll|. It has been shown that there is no Hawking-Page phase transition for the 
Schwarzschild-adS black hole whose horizons have vanishing or negative constant curvature |12l | . 

Thermodynamic properties of black holes in GB gravity with NLED have been studied before [l^. In this paper, 
we will obtain the black hole solutions of GB-Maxwell corrected (GB-MG) gravity in arbitrary dimensions and discuss 
thermal instability conditions. 

One may ask the motivation for considering the GB-MC gravity. As we know, the Maxwell theory has acceptable 
consequences, to a large extent, in various physical areas. So, in transition from the Maxwell theory to NLED, it is 
allowable to consider the effects of small nonlinearity variations, not strong effects. In other words, in order to obtain 
physical results with experimental agreements, one should regard the nonlinearity as a correction to the Maxwell 
field. Eventually, we should note that, although various theories of NLED have been created with different primitive 
motivations, only their weak nonlinearities have physical and experimental importance. So the effects of nonlinearity 
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should be considered as a perturbation to Maxwell theory. In addition, in a gravitational framework the GB gravity 
is a natural generalization of Einstein gravity (not a perturbation in general) in higher dimensions. One may regard 
GB and MC as the corrections of an Einstein-Maxwell black hole. 

The layout of the paper will be in this order: In Sec. m we will introduce the structure of action which contains 
GB gravity coupled with MG. Then, we will solve the field equations and discuss the geometric properties of the 
topological black holes. We will study thermodynamic properties and stability of the asymptotically flat spacetime 
in Sec. mi In the next section, we will generalize our solutions to rotating black branes with adS asymptote and 
calculate thermodynamic and conserved quantities of rotating solutions. Last part is devoted to investigating the 
stability of the solutions in canonical and grand canonical ensembles. Conclusions are drawn in the last section. 


II. TOPOLOGICAL BLACK HOLES 

The Lagrangian of Einstein-GB gravity coupled to NEED can be written as 

Ltot = Le — 2A +aLcB + L{F), (1) 

where the Lagrangian of Einstein gravity is the Ricci scalar, Le = R, and A is the cosmological constant. In the third 
term of Eq. m, a is the GB coefficient with dimension (length)^ and Lcb is the Lagrangian of GB gravity, 

Lgb = RabcdR^'^^^ - ^RabR^'^ + R^. ( 2 ) 

The last term in Eq. ([1]) is the Lagrangian of NLED, which we choose a quadratic correction in addition to the 
Maxwell Lagrangian 

L{F) = -F + PR^ + 0{I3'^), (3) 

where (3 is the nonlinearity parameter and the Maxwell invariant F = FabF°'’’, where Fab = daA}, — dbAa is the 
electromagnetic field tensor and Aj, is the gauge potential. For the vanishing the nonlinearity parameter, this La¬ 
grangian yields the standard Maxwell theory, as it should. Regarding the gauge-gravity Lagrangian o and using the 
variational method, one can obtain the following field equations: 

Gf, + Agab + = ^gabHR) - ^L^Fa^F^, (4) 

da {V^LeF‘^^) = 0 , ( 5 ) 

where Gf), is the Einstein tensor, = 2 {RacdeR-l'^'' - 2RacbdR'"^ - 2RacRl + RRab) - ^LoBgab and Ljr = . 

Now, we are interested in topological static black hole solutions; therefore, we take into account the following static 
metric: 


ds^ = -f{r)dt^ + + r^dnl, (6) 

f[r) 

where dfl^ represents the line element of an (n—I)-dimensional hypersurface with the constant curvature {n—l){n—2)k 
and volume Vn-i with the following explicit form: 



n —1 2—1 

^^1 + S n k = 1 

i=2 j=l 

n— 1 2—1 

dOf + sinh^ OidO^ -I- sinh^ X) 0 k = —1 . 

i=3 j=2 

n — 1 

X k = o 

2=1 


( 7 ) 


Since we are looking for the black hole solution with a radial electric field, we know that the nonzero components 
of the electromagnetic held are 


Ftr — —Frt- 


( 8 ) 
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One can use Eq. ([5]) to obtain the explicit form of Ftr with the following form: 


Ft, = + 0(/32) 


(9) 


which for small values of /3 reduces to Maxwell electromagnetic field tensor. 

We find that the nonzero independent components of Eq. ([4]) can be written by applying electromagnetic field 
tensor dH]): 


eti = F 1^1 + 
2Ar4 






+ 


2g2 




(n — 1) (n — l)r 2"-6 _ ]^y4n-8 




( 10 ) 


egg = r 


1 + 


2(1-/)a' 


+2r(n — 2) 


/" - 2a'r2/'" 

W-2) 


r 


-a{n - 4)(n - 5)(1 - ff - (n - 2)(n - 3)r^(l - /) 


+2Ar‘^ - 


2q2 


+ 


12g4/3 


0(13^). 


( 11 ) 


where a' = {n — 2)(n — 3)a 

It is straightforward to show that the following metric function satisfies all of the field equations simultaneously: 


with 


'l'(r) = 1 + 


Sa' 


n{n — 1) 


A 


/(r) = fc + ^ (l - , 

n{n — l)m ’ 


nq 


2nq‘^P 


2r" 


(n — 2)r^” ^ ^(3n —4) 


oip^), 


( 12 ) 


(13) 


where m is an integration constant that is related to mass. One can see that for small values of /3, the metric function 
reduces to usual GB-Maxwell gravity. Expansion of the metric function for small values of the GB parameter leads to 


/ (r) = Jem - 


4g^ 


a , (fc - IemY , , 
p H-^- a + 


‘(fc “ Jem) 


(n — 1) (3n — 4) ® r 

where the metric function of Einstein-Maxwell gravity is 

2Ar^ 


fEM = k — 


{n — 1) (3n — 4) r 


2g2 


,4n—+ (14) 


(15) 


n(n — 1 ) 7 -" 2 (^77 _ 2 ) 4 

In order to consider the asymptotic behavior of the solution, we put m = q = 0 where the metric function reduces 


to 


fir) = k + —ll-\h 


8a'A \ 
i{n-l) ) ■ 


(16) 


This result puts a restriction on a' which states that for a' < , the metric function will be real asymptotically. 

On the other hand, for A = 0, asymptotically flat solutions are available only for /c = 1 . 

The next step is to look for the singularities. It is a matter of calculation to show that the Kretschmann scalar of 
metric is 


= /"' + 2(n - 1)^ + 2(n - l)(n - 2)^, 


(17) 
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where its series expansion for small and large values of r will be 

lim oc 

r —^0 


( 18 ) 


lim ^ + 

^—>-oo na a \ 


(n — 1) (n + 2) 


(19) 


n — 2 


8 a'A 


^ = 


- 1 , 


n n(n — 1) 
which, for small values of the GB parameter, it will lead to 

lim ^ 


(n + 2)t 


(n-l) 


A + O (a') 


( 20 ) 


( 21 ) 


Eq. (fT5)) - (l?T1) show that there is an essential singularity located at r = 0 and the asymptotical behavior of the 
solutions are adS with an effective cosmological constant. If this solution contains the horizon then our metric 
function is interpreted as a black hole. In order to investigate the existence of the horizon, one should find the root of 
^rr _ suppose that the metric function is positive for large and small r, by considering the possibility 

of the existence of only one extreme root r_|_ = rgxt, we know that /(r) has a minimum at r = rext- So, we can 
investigate the roots of /', 


k(n — l)(n — 2)r] 


.4n—6 


a'{n — 4) 


+ 1 


ext { 2 

L {n - 2)ri„,t 

where, for example, in the Ricci-flat case (fc = 0), we obtain 

[g2A2"-3(±Vl-f 8A/1 - 1)] 


+ 4g4/3 - 2r2"-2g2 _ 2rt:,-^A + O {(3^) = 0, 


'^ext — 


A 


0{p^). 


It is matter of calculation to show that the extremal mass for fc = 0 is 

2g2 2rS,,A 


'^ext — 


rLf(n-l)(n-2) n(n - 1) r3-"(3n - 4)(n - 1) 


O(^2). 


( 22 ) 


(23) 


(24) 


We should note that the obtained solutions may be interpreted as black holes with inner and outer horizons 
provided m > niext, an extreme black hole for m = niext, and naked singularity otherwise. One can obtain the 
Hawking temperature by surface gravity definition with the following explicit form: 


T = 


k{n - l)(n - 2)rt"-'^ (l + - 2rt"-^A - + 4q^P 


47 r(n - (l + 

It is worthwhile to mention that extremal black holes have zero temperature. 


O (P^). 


(25) 


III. THERMODYNAMICS OF ASYMPTOTICALLY FLAT BLACK HOLES (A: = 1 & A = 0) 


In this section we are interested in the thermodynamics of asymptotically flat solutions. Using the series expansion 
of the metric function m with k = 1 for large values of distance, we obtain 


f{r) = 1 - 


m 

j^n —2 


{n 


2g" 

l)(n — 2)r2"“4 


+ 0 



(26) 


which shows that the solutions are asymptotically flat. 
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According to area law, the entropy of black holes is one-quarter of the horizon area. This relation is acceptable for 
Einstein gravity, whereas we are not allowed to use it for higher derivative gravity. For our GB case we can use the 
Wald formula for calculating the entropy of asymptotically flat black hole solutions 

S' = i y (l 2altj (27) 

where R is the Ricci scalar for the induced metric "fab on the (n — 1) dimensional boundary. Calculations show that 
one can obtain 


S = 


Vn-l 

4 



2 (n — 1) 

(n - 3)4 ) 


^n — 1 


(28) 


which confirms that asymptotically flat black holes violate the area law. 

Considering the flux of the electric field at infinity, one can find the electric charge of black holes with the following 
form: 


Q = 


_q_ 

47r 


(29) 


Next, for the electric potential, $, we use the following definition 




q 

(n-2)rr^ 


4(n- 2)g4 \ 

(3n-4) rf-^J 


+ 0(4), 


(30) 


where is the null generator of the horizon. It is notable that although the electric potential depends on the 
nonlinearity of electrodynamics, the electric charge does not. 

In order to calculate the finite mass of the black hole we use the ADM (Arnowitt-Deser-Misner) approach for large 
values of r, which will result in [l^ 


M = ^ m (n — 1). (31) 

IGtt 

We should note that one can obtain m from /(r = r+) = 0, so the total mass depends on GB and NEED parameters. 
Next, by using Eqs. (HH), (Ell), EH) and considering M as a function of the extensive parameters S and Q, we have 


M{S,Q) 


IGtt 


(n 


1 ) 



327r2Q2 10247r4Q4 

(n - 2)4"-^ “ (3n - 4)4"“® 


(32) 


Now, we calculate the temperature and electric potential as the intensive parameter with the following relation 


T = 





(33) 


which are the same as the ones that were calculated in Eqs. EH and EHl). Thus, the conserved and thermodynamic 
quantities satisfy the first law of thermodynamics, 


dM = TdS + 


(34) 


A. Stability of the solutions 

Here, we investigate the thermodynamic stability of charged black hole solutions of GB gravity with NEED. In 
the grand canonical ensemble, the thermodynamic stability may be carried out by calculating the determinant of the 

Hessian matrix of M{S, Q) with respect to its extensive variables W, = (^ gx ■ In the canonical ensemble, 

the electric charge is a fixed parameter, and therefore the positivity of the heat capacity Cq =T+I is sufficient 
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to ensure the local stability. Since the physical black hole solutions have positive temperature, it is sufficient to check 

the positivity of (^) =(|^) /(|£) 


/d^M\ 


An-8 


(n — 1) {n — 2 ) 


2 (n —4)a'^ (n — 8 )a'r^ 4 

(n- 2 ) +"+. 


- 2 r ^”"\2 [2a'(2n - 5) + (2n - 3)r^] + 2Pq^ [Aa 


(4n - 7) + ( 8 n - 10)r2 ] } / tt (n - 1)^ (2a'+ + O (/S^) . 

One may study the behavior of for small values of the GB parameter with the following form: 


(35) 


d^M\ (n-l)(n-2)4”“'‘-2(2n-3)r^"g2 4 (4n - 5) 

~^)q ~ (n - 1)" {n-l fr0-^n^ 

, [(5n-4)r4"-16r2"+V] , , 64g4 

(n - 1) “ (n - 1) 

+O{a'^0). 


In addition, it is worthwhile to mention that for small and large values of r_|_, we obtain 


/92m\ 


V^y, 


= 


4 ( 4 n— 5 )/ 9 q' 
7r(n—1)^ 


^ < 0 , 


Small r-j- 


Large r+ 


27r(n—® 

2(2n-3)q^ ^ ri 

,y(„_l)2^a"-4 > U, 


of r+. 

we obtain 

, a' 


o,p 


0 

a' 

= 

o,p 


0 

a' 


o' 

= 

0 ’ 

a' 

= 

o,p 

= 

0 


n — 2 


r(n — l)r 


< 0 , 


(36) 


(37) 


(38) 


which indicate that for nonzero a' and /3, asymptotically flat black holes with a large or small horizon radius are not 
stable. In other words, there is an upper limit, r+niax, as well as a lower limit, r+min, for the asymptotically flat 
stable black holes (r+min < r+ < r+max)- Although this result does not depend on the value of the GB parameter, 
asymptotically flat black holes with a small horizon radius may be stable in the presence of a pure Maxwell field 
(absence of nonlinearity of electrodynamic correction). In other words, taking into account thermal stability in the 
canonical ensemble, the nonlinearity of electrodynamics has a reasonable influence on a small horizon radius. 

In order to investigate the effects of GB gravity and NLED on the local stability for arbitrary r+, we plot Figs. [T] 
-E Here, we consider left figures which are appropriate for the canonical ensemble analysis. 

In Figs. [Hand[2] we investigate the effects of considering GB gravity. One can find that when the GB parameter is 
less than a critical value (a' < a),), r+ has a real positive root at rg. In this case, (^^0^ > 0 for ro < r+ < r+max- 

In other words, small black holes are not physical and large black holes are not stable. Moreover, we find that 
ro(r+„iax) decrease (increase) when we increase a'. In Fig. |2l we consider a' > a(,, in which T+ is a positive definite 
function of r^. In this case one can obtain stable solutions for r+min < < r+max, in which confirms that small 

and large black holes are unstable. We should note that increasing a' leads to increasing (decreasing) the value of 
?'+min (?'+max) and, therefore, decreases the range of stability. 

Besides, we plot Figs. |3]and|l]to analyze the effects of NLFD. We find that for the fixed values of q, n and a', there 
is a critical value for the nonlinearity parameter, Pc, in which for /3 < pc, the temperature is positive for > tq. 
This case is the same as that in Fig. [I]with a' < a'^ and we may obtain asymptotically flat stable black holes for 
To < r+ < r+max- In Fig. |3]we set P > Pc to investigate the stability condition for positive definite temperature. 
One finds there are lower and upper limits for the horizon radius of stable black holes. It means that for P > Pc, 
asymptotically flat black holes are stable when r+min < r+ < r+max- 

We study the effects of electric charge in Figs. |5]and|ni These figures show that for fixed values of a' and p, there is 
a critical value for the electric charge, qc, in which the temperature is positive definite for q < qc- When q > qc, there 
is a real positive root (rg) for the temperature which is an increasing function of q. We should note that the general 
behaviors of Figs. [5] and [6] are, respectively, the same as those in Figs. [T]and|2l In other words, asymptotically flat 
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FIG. 1: Asymptotically flat solutions: (left figure) and |H^q| (right figure) versus r+ for n = 5, g — 1, fl = 0.001 

and a = 0.1 (solid line), a = 0.3 (dotted fine) and a = 0.5 (dashed line), ’’bold lines represent the temperature” 



( 2 \ II 

^0- ] (left figure) and 1—(right figure) versus r+ for n = 5, 5 = 1, /3 = 0.001 and 
a = 0.6 (solid line), a = 0.8 (dotted line) and a = 0.1 (dashed line), ’’bold lines represent the temperature” 


black holes are stable for r+min < r_|_ < r+max and one should replace r+min with tq for q > Qc- It is worthwhile to 
mention that both r+min and r+max are increasing functions of q. 

Finally, considering Figs. [T]-[ni we should note that, in canonical ensemble, asymptotically flat black holes are 
stable for r+min < r+ < r+max) in which one must replace r+min with tq when T+ has a real positive root {q > qc or 
13 < Pc or a' < a'P). 

Next, in the grand-canonical ensemble, we calculate the determinant of the Hessian matrix of the asymptotically 
flat black holes. After some algebraic manipulation, one obtains 

|H^q| = 4 {2q\in - 4)rt"-" [r^ - 2a'] + 

4g'^r^"-2 [2a'(lln - 24) - (7n - 8)r2 ] 































FIG. 3: Asymptotically flat solutions: (left figure) and IH^^qI (right figure) versus r+ for n = 5, g = 1, a = 0.1 and 

P = 0.001 (solid line), p ~ 0.003 (dotted line) and p = 0.005 (dashed line), ’’bold lines represent the temperature” 



FIG. 4: Asymptotically flat solutions: ^j (left figure) and |H^q| (right figure) versus r+ for u = 5, g = 1, a = 0.1 and 
p = 0.01 (solid line), P = 0.03 (dotted line) and p = 0.05 (dashed line), ’’bold lines represent the temperature” 



(39) 







FIG. 5: Asymptotically flat solutions: ^ j (left figure) and |H^q| (right figure) versus r+ for n = 5, /3 = 0.001, a 
and q = 1.2 (solid line), q — 1.3 (dotted line) and q — 1.4 (dashed line), ’’bold lines represent the temperature” 




FIG. 6: Asymptotically flat solutions: ^ ^ (left figure) and |H^q| (right figure) versus r+ for n = 5, (3 = 0.001, a 

and q = 0.5 (solid line), q — 0.7 (dotted line) and q — 0.9 (dashed line), ’’bold lines represent the temperature” 


In addition, it is worthwhile to mention that for small and large values of r_|_, we obtain 


ttM 

Small 


< 


4{lln-24)/3q‘^ . „ 

7r(n-l)2(3n-4)a'2r“"-"’’ 

16(7n-8)/3g^ I 

r“"-“(n-l)2(n-2)(3n-4) 

_^ Q 

_^ n 

(n-l)2(n-2)r^"-® ’ 

4 


0 , /3 0 

a' = 0,/3 ^ 0 
a' ^ 0,/3 = 0 


a' = 0,/3 = 0 


7r(n — l)r” 


jjM 

Large 


< 0 , 
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which indicate that, for nonvanishing a' and /3, asymptotically flat black holes with a large horizon radius are not 
stable. In other words, there is an upper limit, r+maxi for the asymptotically flat stable black holes (r_|_ < r+max)- 
In the absence of the GB gravity case (o' = 0,/3 ^ 0), there is a lower limit for stable solutions. It means that for 
a' — Q the results of the stability conditions for the canonical and grand canonical ensembles are identical and both 
ensembles show that small and large black holes are unstable. We should note that in the presence of GB gravity 
thermal stability is ensemble dependent. These results are shown in Figs. [T]"|6l 

Now, we need to discuss ensemble dependency 1^. As we know in the canonical ensemble the internal energy is 
allowed to fluctuate with fixed electric charge and, therefore, the black hole and the heat reservoir remain in thermal 
equilibrium with a certain temperature. While in the grand canonical ensemble the black hole is in both thermal and 
electrical equilibrium, with its reservoir held at a constant temperature and a constant potential. In other words, 
different boundary conditions lead to different ensembles. 

In the usual discussions of the stability criterion of black holes, one expects ensemble independence of the system. 
Indeed, different ensembles which imply different boundary conditions, should lead to similar stability conditions in 
the usual thermodynamical systems. 

Ensemble dependency of a system may come from two subjects. One of them is real ensemble dependency, which 
may occasionally occur in some thermodynamical systems. The existence of ensemble dependency was seen in the usual 
thermodynamical systems 0- Another one, which is not real, can be removed by improving our thermodynamic 
view point. In other words, our lack of knowledge may lead to ensemble dependency and we should improve our 
thermodynamical understanding to obtain ensemble independency. In our black hole model (with a spherical horizon), 
we find that the nonlinearity of electrodynamics results in the same changes for the behavior of both canonical and 
grand canonical ensembles. One can see the ensemble dependency comes from the contribution of Gauss-Bonnet 
gravity. This means that in the absence of the GB parameter, both ensembles have the same stability conditions. This 
shows that the GB parameter makes more of a contribution to thermodynamical behavior of the black hole systems. 
In other words, one may consider the GB parameter not only as a fixed parameter, but also as a thermodynamical 
parameter 0. By doing so, the Hessian matrix for this system will be modified and, therefore, we may expect to 
see that this modification solves the ensemble dependency of the black hole system. 

Another way to solve ensemble dependency is through considering the fact that thermodynamical systems described 
by a thermodynamical potential must be invariant under the Legendre transformation. To do so, one can use the 
method that was introduced in and can use geometrothermodynamics to solve ensemble dependency of the 
system p^[l^. 


IV. ASYMPTOTICALLY ADS SPINNING BLACK BRANDS (fc = 0 & A / 0) 

In this section, we will investigate rotating adS space-time. In order to investigate the asymptotic behavior of the 
solutions for fc = 0 and A 0, one can use the series expansion of the metric function for a large value of r. We 
obtain 


\ 2Ae£f 2 

j{r) =-^- —r - 


m 


2q^ 


n(n-l) (i + ^yn -2 („_!)(„_ 2) (l + ^y 


O 


.2n-4 


1 


r>2n—2 I ’ 


(42) 


where 


Aeff = 


i(n- 1) 


1 I Sa'A 1 ' 
^ + n(n-l) “ ^ ; 


4a' 


(43) 


Eq. (l42l) shows that for fc = 0 the solutions are asymptotically adS with an effective cosmological constant Aeff. When 
constructing a rotating space-time one can apply the following transformation in the static space-time with fc = 0 


t 


it 

^ t. 


(44) 


Using this transformation, the metric of {n-\- l)-dimensional asymptotically adS rotating space-time with p rotation 
parameters is 


ds^ 


-f{r) 


P 

'E.dt — aid(l)i 


i=l 


dr^ 

f(r) 


2 P 

^ ^ {a^dt - ^Pd(j)i)^ 


2 
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^ iaid(j)^ - ajd4>if + r^dX'^, 

i<j 


(45) 


where .n = yl + X)r=i W Euclidean metric on the (n — 1 — p)-dimensional submanifold. The forth 

term in Eq. (1451) comes from the fact that for rotating space-time with more than one rotating parameter, we should 
consider cross terms associated with rotating coordinates. The rotation group in {n 1) dimensions is SO{n) and, 
therefore, p < [n/2]. Considering the aforementioned transformation, for the gauge potential, we should write 

= h(r) (S5° — (no sum on i). (46) 


Calculations show that metric function (ED with k = 0 and h{r) = f Ftrdr (using the Ftr calculated in Eq. (0) 
satisfies all of the field equations for the aforementioned rotating adS spacetime. Straightforward calculations confirm 
that there is a curvature singularity located at r = 0 which is covered with the horizon(s). 

Using the analytic continuation of the metric by setting t —>■ fr and > iai, and regularity at the event horizon, 
r+, helps us to obtain the Hawking temperature and the angular velocities of the black branes: 


/'(r-+) 

47r.^ 


An—4: ! 


A + - 2q^l3 

27rS(n - l)r^"-5 


+ o{p% 


(47) 


= (48) 

Eq. (ITTl) shows that, unlike black hole solutions with a spherical horizon, the temperature of the aforementioned black 
brane with a flat horizon does not depend on GB gravity. 

Next, we calculate the electric charge and potential of the solutions. The electric charge per unit volume 14,_i can 
be found by calculating the flux of the electric field at infinity, yielding 


Q = 


47r 


(49) 


On the other hand, because of the applied transformation and changes in metric we now have Killing vectors in the 
form ot X = dt + X)r=i which is the null generator of the horizon. The electric potential is obtained as 




g / _ 4(n-2)g^/3 \ 
S(n - 2 )r !(-2 ( 3 n- 4 )r 2 "- 2 y 


+ 0(/3")- 


(50) 


Now, we desire to calculate the entropy, the angular momentum and the finite mass to check the first law of 
thermodynamics. In general, the action and the conserved quantities of the space-time diverge when evaluated on the 
solutions. In order to overcome this problem and due to the fact that our space-time is asymptotically adS, we can 
use the counterterm method to calculate the finite action and the conserved quantities. One can show that for the 
obtained solutions with flat boundary, Rabcd{l) = 0 , the finite action is 


.Infinite — F Ih F let ; 


(51) 


where 


Ig = - 


IGtt , 
1 




'M 




let = - 


h 




'dM 


n — 1 

^eff 


(52) 

(53) 

(54) 


where 7 and K are, respectively, the trace of the induced metric, Xab and the extrinsic curvature, K°‘^, on the boundary 
dM, Gab is the Einstein tensor calculated on the boundary, J is the trace of 


Jab = liKedK^’^Kab F 2KKaeK^ - 2KaeK^'^Kdb - K^Kab) 


(55) 
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and /eff is a scale length factor that depends on I and a, which must reduce to Z as a vanishes. It is worthwhile to 
mention that, for a spacetime with zero curvature boundary, let has exactly the same value as that of the Einstein 
gravity in which I is replaced by leg. Using the Brown-York method of a quasilocal definition with Eq. (l5T|) - (l54)) . one 
can introduce a divergence-free stress-energy tensor as follows 


T 


ab 


1 

Stt 



-h 2a(3 




leS 


(56) 


Then the quasilocal conserved quantities associated with the stress tensors of Eq. dSSl) can be written as 


QiO = J d”-VV77’abn“e', (57) 

B 

where n“ is the timelike unit normal vector to the boundary B. Taking into account Eq. (1571) with ^ = d/dt as a 
Killing vector, one can calculate the mass per unit volume Ui-i as 

M = (nS^ - l) , (58) 

where 

m = m{r = r+). (59) 

Eqs. (|5^ and (|5^ indicate that, unlike the asymptotic flat solutions with a spherical boundary, the finite mass does 
not depend on the GB parameter for the boundary flat rotating solutions. 

Considering another Killing vector C = d/d(j)i of rotating spacetime which is related to angular momentum, one 
can calculate the angular momentum per unit volume, 

Ji = -^nmEtti. (60) 

IQtt 

We should note that for static solutions (a^ = 0), the angular momentum vanishes, which confirms that the a^’s are 
the rotational parameters of the space-time. 

As we mentioned before, the Wald formula can be applied to asymptotically flat space-time. Here, we encounter 
asymptotically adS solutions, so we calculate the entropy through the use of the Gibbs-Duhem relation. 


S = - ) - -^finite- 


(61) 


First we calculate the finite action /finite for the rotating metric. We find that the finite action per unit volume may 
be written as 


/finite — 


87rn(n — 1)T4. 


-Ar" + 


nq 


2nq^f3 


(n — 2)r 


(3n — 4)r 


,3n—4 
+ 


om 


Using the finite conserved and thermodynamic quantities with the finite action, we obtain 


(62) 



(63) 


which confirms that the entropy obeys the area law for asymptotically adS black branes with zero curvature horizon. 

Now we want to check the first law of thermodynamics. To do so, we obtain the mass as a function of the extensive 
quantities S', J and Q. Using the expression for the electric charge, the mass, the angular momentum and the entropy 
given, respectively, in Eqs. dSl), (EHD, dini) and (1551) and the fact that /(r = r+) = 0, one can obtain a Smarr-type 
formula as 


M(S, /,Q) 


[nZ - 1) J 
nl^jZ{Z - if 


(64) 
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where J = |J| = \/Yli=i Z = is the positive real root of the following equation 

(3n — 4) {3n — 4) AS"* \ 


2^nl^/Z{Z -1) (^7r'*Q'*/3 - 


2(n-2) 


2n J 


+ 




S'^Ztt (n — 1) (3n — 4) J = 0. 


(65) 


Taking into account S, J and Q as the extensive parameters of M, we can define the intensive parameters conjugate 
to them as 


/aM\ _fdM'' 


( 66 ) 


It is a matter of calculation to show that the intensive quantities calculated by Eq. (IS51) coincide with Eqs. (I47L 
and dsni). Therefore, the first law of thermodynamics is satisfied 


dM = TdS + + <^dQ. 


(67) 


A. Stability of the solutions 


Now, we are in a position to calculate the heat capacity and the Hessian matrix to check the local stability of these 
solutions in context of canonical and grand canonical ensembles. For canonical ensembles, where electric charge and 
angular momenta are fixed parameters, the positivity of ^ is sufficient to ensure the local stability. Therefore, 

we obtain 


1 ^ 


2 r;"Ai 


4(n — 2)^(7*r|‘*A2/3 


d^M 

dS^ J jQ (’^ ~ 1) l/’O' _ X) (3n — 4) cr'0^ 


O (/?'), 


( 68 ) 


Ai = 


(3cr - n) 2 (3tT - n^) Aq"^ [{n + 2)cr - n^] 


3n—6 


+ 


1—4 


+ 


(69) 


A 2 = 


[{5n^ - 42n + 40)S2 + {7n^ + lln - 20)] A^ 


^n — 6 


2n [(13^2 - 50n + 40)S2 - (n^ - 19n + 20)] Aq^ 


{n — 2)r 


3n—8 
+ 


+ 


[(17n — 20 )ct — n(5n — 6)] q^ 


(n-2) 


2„5«-10 


(70) 


ij} = \jiq^ - A{n - 2)r^" ^] , 
a = [{n- 2)S^ + 1 ] . 


(71) 

(72) 


As we mentioned before, in the grand canonical ensemble, the positivity of the determinant of the Hessian matrix 
of M{S, Q, J) with respect to its extensive variables Xi, = (^ gx gx ) ’ sufficient to ensure the local stability. 

It is a matter of calculation to show that the determinant of H^g j is 

A,Q,j\ — ^(9 — Ar_|_ j — 

A^iq'^ [3(n — 2)^A^r*"“‘* + 3n(n — l)q^ — 2(n — 2)(3n — 2)Ag^r^"“^] /3 


(3n - 4) [V - (n - 2) Ar^"-^] 


0(/3'), 


(73) 
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where 

647rrt”^" [nq'^ - {n - 2) Ar^"-^] 
[(n-2)S2 + l] S6Z2 


(74) 


Regarding Eq. 


(H71) with the mentioned 



and also 


IwM I 


we find that, unlike the asymptotically flat 


case with the spherical horizon, neither the heat capacity nor the determinant of the Hessian matrix depend on the 
GB parameter for asymptotically adS rotating solutions with zero curvature horizon. Therefore, we expect to obtain 
the same results in both canonical and grand canonical ensembles. Following the method of the previous section and 
regardless of the values of n, q, A, S and a, one finds 


and 


d^M\ 


5S'2 ^ 

J-.Q 

Small 




dS^ ^ 

J,Q 

Large 


tjM I 

S,Q,J \ Small 


1 

1 Large r+ 


4[(5n^-23ra+20)(H=-l) + (12n=*-31n+20)E:^](3"‘/3 a / n 

7rH2(„-l)(3n-4)<Tr5"-“ ^ 

2{3[(»-2)5^ + i]-»K . n 

7rH2(„-l)[(n-2)H2 + l]r^"-* ^ ^ 

2 [(n^ — 4) (S^ — l) + n — 2] A ^ 

' 7rS2 {n - 1) (n - 2) [(n - 2)S2 + 1] r”+2 ^ 


_ 7687r(n-l)g^/3 ^ R C) 

6477 


n[(ri-2)H2 + l]H6i2ri)" 

647r 


(n - 2)Z2S6 [{n - 2) S2 + 1] r* 


— >0, ^ = 0 ’ 

> 0 . 


3n-4 


(75) 

(76) 


(77) 

(78) 


Both ensembles confirm that, in the presence of NEED (/3 ^ 0), although the black branes with small r+ are unstable, 
large black branes are stable. It is notable that the instability of the small black branes is due to the presence of the 
NEED and, in the absence of the nonlinearity effect, small black branes are stable. 


V. CLOSING REMARKS 

In this paper, we regarded both the gravity and electrodynamic string corrections of Einstein-Maxwell gravity to 
obtain black hole solutions with spherical, hyperbolic, and fiat horizon topology. 

At the first step, we focused on asymptotically flat solutions. We used the Wald formula, the Gauss law and the 
ADM approach to calculate entropy, electric charge and finite mass, respectively. We checked that the conserved 
and thermodynamic quantities satisfied the first law of thermodynamics. Then we investigated the thermodynamic 
stability of the solutions in both canonical and grand canonical ensembles. Taking into account the canonical ensemble, 
we found that for nonzero a' and /3, asymptotically flat black holes with a large or small horizon radius are unstable. 
This means that, in canonical ensemble, asymptotically flat black holes are stable for r+min < r+ < r+max, in which 
one must replace r+min with ro (the largest root of T+) when T+ has a real positive root (q> qcOT j3 < /3c or a' < a(,). 
Moreover, we found that different values of a, /3 and q can change the values of r+min and r+max- 

Then, we investigated the stability conditions in the grand canonical ensemble. We showed that there is an upper 
limit, r+max, for the asymptotically flat stable black holes (r+ < r+max)- We found that, although for a' = 0 the 
results of the stability conditions for canonical and grand canonical ensembles are identical, these ensembles have 
different consequences in the presence of GB gravity. In other words, we noted that in the presence of GB gravity, 
thermal stability is ensemble dependent. 

In the next section, we considered the Ricci flat solutions with an adS asymptote and produced a rotating spacetime 
by using an improper local transformation. In addition, we calculated the conserved and thermodynamic quantities 
for asymptotically adS black branes which satisfy the first law of thermodynamics. Gonsidering the thermodynamic 
instability criterion in canonical and grand canonical ensembles, we found that neither the heat capacity nor the 
Hessian matrix depend on the GB parameter. Therefore, both ensembles have identical stability conditions. We 
found that, although the black branes with small r+ are unstable, large black branes are stable (unlike asymptotically 
fiat static large black holes which are unstable). It is notable that the instability of small black branes (holes) is due 
to the presence of NEED and, in the absence of the nonlinearity effect, small black branes (holes) are stable. 

To conclude, we found that for the Ricci-flat solutions, we obtained the same conditions for stable black holes in 
both canonical and grand canonical ensembles. In other words, in this case, we took into account S, Q and J as 
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the thermodynamical extensive parameters, correctly, and obtained the same results for both ensembles. It is easy 
to show that, for rotating Ricci-flat solutions, one may obtain ensemble dependency if one, imprecisely, considers S 
and Q as the set of extensive parameters (regards J as a dynamical parameter, not a thermodynamic one). As one 
can confirm, unlike Ricci-flat solutions, the GB parameter in the spherical horizon of GB black holes contributes to 
finite mass, temperature and, consequently, heat capacity. In other words, the GB parameter contributions lead to 
ensemble dependency in spherical horizon black holes. This means that a GB generalization of Einstein gravity not 
only affects gravitational properties, but also thermodynamic aspects of the spherically symmetric black holes. 

Finally, we should note that the modifications of the thermodynamic instability criterion in the presence of the GB 
gravity depend on the choice of the ensemble. One may consider the GB parameter as a thermodynamic variable to 
remove the ensemble dependency. In addition, it is worthwhile to mention that we can regard asymptotically adS 
black holes with spherical topology to investigate P — V criticality in the extended phase space of the solutions by 
calculating the Gibbs free energy. These works are under examination. 
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